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ABSTRACT 


We compare the dynamics of satellite galaxies in the Sloan Digital Sky Survey to 
simple models in order to test the hypothesis that a large fraction of satellites co-rotate 
in coherent planes. We confirm the previously-reported excess of co-rotating satellite 
pairs located near diametric opposition with respect to their host, but show that this 
signal is unlikely to be due to rotating discs (or planes) of satellites. In particular, no 
overabundance of co-rotating satellites pairs is observed within ^ 20° — 50° of direct 
opposition, as would be expected for planar distributions inclined relative to the line- 
of-sight. Instead, the excess co-rotation for satellite pairs within ^ 10° of opposition is 
consistent with random noise associated with undersampling of an underlying isotropic 
velocity distribution. Based upon the observed dynamics of the luminous satellite 
population, we conclude that at most 10% of isolated hosts harbor co-rotating satellite 
planes (as traced by bright satellites). 

Key words: Local Group - galaxies: formation - galaxies: evolution - galaxies: dwarf 
- galaxies: star formation 


1 INTRODUCTION 

Within the ACDM paradigm, the growth of cosmic struc¬ 
ture proceeds as overdensities collapse into dark matter ha¬ 
los, which eventually serve as the sites for galaxy formation 


(White & Rees 1978 Blumenthal et al. 1984 Davis et al. 


1985). Over time, the hierarchical accretion and merging of 


halos drives the development of substructure, such that some 


halos reside within the bounds of larger parent halos (Moore 
et al.|1998 ). The galaxies hosted by these subhalos are typ¬ 


ically referred to as satellites and are important probes of 
the evolution of substructure, as they serve as tracers of dark 
matter on small scales. 

High-resolution A-Body and hydrodynamic simulations 
confirm this picture of hierarchical structure formation, 
while also making predictions regarding the properties of 
subhalos and the satellite galaxies they host. In particular, 
simulated subhalos are not isotropically distributed with re¬ 
spect to their parent dark matter halo. Instead, simulations 
across a broad range of mass scales predict that satellite 
galaxies should preferentially lie along orbits aligned with 


the major axis of the host halo (e.g. van den Bosch et al. 
|1999| |Knebe et al.||2004l [Libeskind et al.||2005| [Kang et al. 


* e-mail: johnip@uci.edu 
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2007 Lovell et al.|2011 Wang, Frenk Cooper|20l3 ). Two 


possible physical drivers are often associated with this pre¬ 
dicted alignment of substructure with the shape of the larger 
gravitational potential: [z] preferential destruction (or sup¬ 
pression) of satellites on orbits anti-aligned with the halo’s 


major axis (Zaritsky & Gonzalez 

[1999 

Peharrubia, Kroupa 

& Boily||2002 Pawlowski et al.| 

12012 

) or [a] accretion of 


satellites along preferred directions, perhaps associated with 


2011 ). 


large-scale filaments (Zentner et al. 2005 Libeskind et al 


Observations of galaxies in nearby groups and clusters 
largely support the predicted anisotropies found in simu¬ 
lations, such that satellites in massive dark matter halos 
are preferentially aligned with the major axis of the cen¬ 
tral galaxy and with the larger-scale, filamentary structure 


(e.g. West & Blakeslee|2000 Plionis et al.|2003 

Faltenbacher 

|et al. 2007 Hao et al. 2011 Tempel et al. 

2015). When 

pushing to lower-mass, more-isolated halos, 
on large spectroscopic samples similarly find 
preferentially reside along the major axis of 
type) hosts, while the distribution of satellite! 
(or late-type) hosts is consistent with being iso 

studies based 
that satellites 
red (or early- 
3 around blue 

tropic (|Brain- 

erd|2005 |Sales & Lambas|2004 |2009[|Yang et al.|20061 |Az- 

zaro et al.|2007||Bailin et al.|2008| but see also|Zaritsky et al. 


1997). This apparent lack of spatial anisotropy for satellites 


of late-type hosts is potentially driven by random misalign- 
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2 Phillips et al. 


merit between the major axis of the host’s disc and the dark 
matter halo, such that the satellites may be aligned with the 
latter but not the former ( Bailin et al.|2005 Libeskind et ah 


2007 Deason et al.||2011 ). 

In contrast to the satellites of comparable star-forming 
hosts in the local Universe, observations of the Local Group 
suggest that the spatial distribution of satellites around both 
the Milky Way and M31 are significantly anisotropic. In par¬ 
ticular, the satellites of the Milky Way preferentially reside 
near the northern and southern Galactic poles (i.e. along the 
minor axis of the Milky Way disc, Holmberg|1969|), possibly 
following a planar arrangement ( Lynden-Bell||197^ [Kunkel 


Demers||1976| |Metz, Kroupa Jerjen||2007| [Pawlowskp 


Pflamm-Altenburg & Kroupa 2012). The satellites of M31 


are similarly anisotropic in their distribution, with a large 


subset belonging to a thin disc or plane ( 

Karachentsev| 1996 

Koch & Grebel||2006| 

McGonnachie & Irwin||2006| |Metz,| 

Kroupa & Jerjen||2009 

Gonn et al.||2013 ). 


When including velocity information, the anisotropy of 
the Local Group satellite distribution becomes even more 
pronounced, with many of the Milky Way satellites follow¬ 
ing polar obits, consistent with a vast, coherently-rotating 
plane ( Metz, Kroupa Libeskind|2Q08| Pawlowski, Kroupa 


Jerjen||2Q13| Pawlowski Kroiipa||2013 ). For M31, a yet 

more-striking planar structure is observed, such that a large 
number of satellites exhibit coherent rotation along the line- 
of-sight to the Milky Way, forming a vast plane with a di¬ 
ameter of ~ 400 kpc and a thickness of less than ~ 15 kpc 
( Ibata et al.|2013 ). While simulations predict that satellites 
should preferentially align with the major axis of the host 
dark matter halo, the strong anisotropies observed for the 
Local Group satellites (especially those around M31) are in¬ 
consistent with the expectations of simulated subhalo pop¬ 
ulations ( Kroupa, Theis fc Boily||2005 Kroupa et al.||20l'0 


Pawlowski et al.||2012[ 2014f ^ata et al.|2014c[ but see also 

Buck, Macci’o Duttonl|2015| who argue that co-rotating 
planar arrangements of satellites are predicted by AGDM). 

The striking nature of the M31 satellite disc has served 
as fuel for many recent studies investigating the possibility 
of similar, strongly-anisotropic satellite distributions around 
galaxies outside of the Local Group, such as the discovery 
of possible planar structure in the satellite distribution of 
the Gentaurus A group ( |Tully et al.||2015' Libeskind et al. 


2015). In particular, recent analysis of satellite pairs in the 
Sloan Digital Sky Survey (SPSS, [York et 2000) points 


towards the possibility of co-rotating planar satellite struc¬ 


tures around nearby massive galaxies (Ibata et al. 


2014a 


hereafter 114); for 20 out of 22 systems, with satellite pairs 
located on diametrically-opposed sides of the host galaxy, 
114 detect co-rotation along the line-of-sight, suggesting that 
thin satellite planes - similar to that of M31 - may be rela¬ 
tively common. Specifically, this result, which is supported 
by an analysis of the spatial positions of photometrically- 
selected satellite samples, indicates that > 50% of the satel¬ 


lite population may reside in thin co-rotating planes (Ibata 


|et al.|[?014b[ although [Gautun et al.|[2015| argue that the 
evidence for the ubiquity of such planar structures is not ro¬ 
bust). Given the scarcity of such structures in modern sim- 
ula tions ([Ibata et al. 


2014c 


Pawlowski & McGaugh 


2014 


see [Gautun et al.||2015| for m argument that the diversity 
of properties of these structures accounts for their perceived 


© . 



Figure 1. Definition of the opening angle (a) between satellite 
pairs (here, SI and S2), as measured with respect to the host 
galaxy. Each satellite pair has a uniquely defined o, ranging from 
0° to 180°, such that satellites on diametrically-opposed sides of 
a host correspond to o = 0°. 


rarity.), the analysis of 114 poses a strong test of the AGDM 
cosmology and thereby warrants further investigation. 

In this paper, we re-examine the kinematic evidence 
for the existence of co-rotating planes of satellites around 
nearby massive hosts by comparing the coherence of line- 
of-sight velocities of observed satellite galaxies to simple 
models of satellite spatial distributions and kinematics. The 
structure of the paper is as follows: in il we discuss the se¬ 
lection of the observational sample and the measured abun¬ 
dance of co-rotating satellite pairs. In E we introduce our 
numerical models and compare the mock observations de¬ 
rived from the models to the observational data. Finally, 
in Q we discuss our results in the context of the search 
for M31-like planes elsewhere in the Universe. Throughout 
our analysis, we employ a A cold dark matter (AGDM) cos¬ 
mology with WMAP7+BAO+i7o parameters Qa = 0.73, 
Urn = 0.27, and h = 0.70 ( Komatsu et al.||2011| ), and un¬ 
less otherwise noted all logarithms are base 10. Throughout 
the paper, we use the terms “[satellite] disc” and “[satel¬ 
lite] plane” interchangeably, referring to co-rotating planar 
satellite configurations. 


2 OBSERVATIONAL DATA 


2.1 Sample Selection 


We draw our observational data from Data Release 7 (DR7, 
Abazajian et al.||2009 ) of the SDSS, making use of the de¬ 
rived data products from the NYU Value-Added Galaxy 
Gatalog (VAGG, Blanton et al.||2005 ) including absolute r- 
band magnitudes (Mr) that are A-corrected to z = 0.1 using 
KGORREGT (Blanton & Roweis 2007). Throughout this 


work, we restrict our analysis to regions of the SDSS where 
the spectroscopic completeness exceeds 70% (i.e. FGOT- 
MAIN > 0.7). We also reject all galaxies with line-of-sight 
velocity errors greater than 25 km/s. 

In selecting our galaxy sample, we adhere closely to 
the procedure of 114. We select a sample of hosts in the 
magnitude range — 23 < Mr < —20 and within a redshift 
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Figure 2. Left: the cumulative number of co-rotating (green line) and counter-rotating (orange line) satellite pairs as a function of 
opening angle (a), including uncertainties based on Poisson statistics. Right: the fraction of co-rotating satellite pairs as a function of 
opening angle, computed in distinct bins of a (bins have width 6°), with error bars derived according to the binomial theorem. There 
is a significant excess of co-rotating satellite pairs over counter-rotating pairs at small a (i.e. for satellites that are nearly diametrically 
opposed to each other). This overabundance of co-rotating pairs suggests the possible presence of coherently-rotating planar satellite 
structures around nearby massive host galaxies (however, see Fig. |^. 


range of 0.002 < z < 0.05. A host is considered isolated if 
there are no brighter objects within 500 kpc (in projection) 
on the sky and within 1500 km/s in velocity space. Only 
isolated systems are retained, reducing the number of hosts 
to 22, 780 isolated galaxies. From this set of host systems, 
we identify galaxies as satellites of a given host if 

(i) their magnitudes fall in the range 

Afr,host + 1 < Mr ,sat < -16, 

(ii) they are located between 20 kpc and 150 kpc from their 

host in projected distance (tiproj), and 

(iii) their velocity offset from the host lies in the range 

25 km/s < \Vsat — Vhostl < 300 km/s x . 

This velocity bound is taken from 114, and is designed to 
reduce the contamination from interlopers in the satellite 
sample. Since our interest is in pairs of satellites, we re¬ 
tain only hosts with two or more satellites. Our final sample 
contains 427 such hosts, with 965 associated satellites. Note 
that individual hosts are allowed to harbor more than two 
satellites; on average, the SDSS hosts (as well as our model 
hosts, see have 2.3 satellites. 

2.2 Co-rotation signal 

In this subsection, we investigate pairs of satellites for evi¬ 
dence of co-rotation with respect to their host. To facilitate 
this, we introduce the parameter o, defined as the angle be¬ 
tween the line extending from one satellite through the host 
and the position vector of the second satellite relative to 
the host, as projected on the sky (see Figure]^. We define 
this opening angle a such that a satellite pair located on 
diametrically-opposed sides of a host will have an opening 


angle of 0°. For the duration of this work, we will refer to a 
satellite pair as “co-rotating” if the satellites have opposite- 
signed (i.e. one + and one -) line-of-sight velocity offsets 
relative to their host and their associated opening angle (a) 
is less than 90°, or if they have same-signed velocity offsets 
relative to their host and their associated a is greater than 
90°. Otherwise, the satellite pair is deemed to be counter¬ 
rotating. 

Figure shows the cumulative number of co-rotating 
and counter-rotating satellite pairs in our sample as a func¬ 
tion of opening angle at a < 35°. At small opening angles 
there is a clear excess of co-rotating pairs, as first reported 
by 114. The overabundance of co-rotating pairs as a func¬ 
tion of a is better illustrated in Figure |^, which shows the 
fraction of satellite pairs that are co-rotating as a function 
of opening angle, computed in distinct bins of a. While the 
surplus of co-rotating pairs at small opening angle is readily 
apparent, at 10° < a < 35° the signal is consistent with 
the sample being divided equally between co-rotating and 
counter-rotating (i.e. a co-rotation fraction of 0.5, as would 
be expected in the absence of any co-rotating structure). At 
first glance, the data would seem to indicate the presence of 
coherently rotating structures that can only be detected at 
small values of a (i.e. co-rotating planes of satellites viewed 
close to edge-on). 

When examining the co-rotating fraction of satellite 
pairs over the full range of opening angles (i.e. 0° < a < 
180°), however, the evidence for planes of satellites is far less 
convincing. In Figure]^ we show (a) the cumulative counts 
of co-rotating and counter-rotating satellite pairs and (b) 
the co-rotating fraction, again in discrete bins of a, over 
the full range of opening angles Q In this light, the excess 

^ Note that the co-rotating fraction as a function of opening angle 
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a (degrees) 



Figure 3. Left: the cumulative number of co-rotating (green line) and counter-rotating (orange line) satellite pairs as a function of 
a, including Poisson errors, over the full 180° domain. Right: the fraction of co-rotating satellite pairs as a function of opening angle, 
computed in distinct bins of a (bins have width 10°), with error bars derived according to the binomial theorem (solid black line). To 
facilitate visualization, a coarser binning is adopted than that employed in Fig.[^; the dash dashed line, however, shows the corresponding 
dependence of co-rotating fraction on a utilizing this narrower binning. Over the full range of opening angles, the observations are largely 
consistent with no excess of co-rotating satellite pairs, suggesting an isotropic velocity distribution. 


of co-rotating pairs at a < 10° seems less likely to result 
from structured, coherent rotation associated with planar 
distributions of satellites. For example, repeatedly resam¬ 
pling 400 satellite pairs placed randomly in phase space 
(as in our “isotropic” model in ^ will frequently produce 
satellite samples with excess co-rotating fractions at random 
opening angles that are, by definition, not indicative of any 
underlying physics. Thus, care must be taken not to over¬ 
interpret the observed overadundance of co-rotating pairs 
at small angles, if indeed it is merely the result of random 
fluctuations associated with undersampling of an underlying 
isotropic distribution. The remainder of this paper will ex¬ 
amine the argument that the excess of co-rotating pairs at 
small a is significant, and indicative of ubiquitous coherent 
co-rotation (similar to that observed for M31’s satellite pop¬ 
ulation) by comparing the SDSS data to statistical models 
of satellite kinematics. 


3 COMPARISON TO TOY MODELS 

In order to gain insight as to whether the data presented in 
Figures and does indeed argue for the existence of co¬ 
herently rotating satellite structures, we compare the SDSS 
data to mock observations of simple, idealized “toy” models 
of satellite systems. These models are not intended to give 


is computed using different binning procedures in Fig. and 
Fig.[^, so as to aid in visualization. Our results do not strongly 
depend on how the data are binned. The increase in co-rotating 
fraction at small opening angle appears less significant in Fig. S’ 
as a result of allowing satellite pairs with opening angles of ~ 10° 
in the innermost bin. 


a detailed description of satellite phase-space distributions 
However, they do provide a meaningful basis for comparison 
to the observational data. We begin by detailing how each 
toy model is constructed: 

(i) Isotropic model — For this model, each host is randomly 
assigned 2 — 5 satellite galaxies]^ The position and velocity 
of each satellite, with respect to the host, is randomly chosen 
to be between 0 and 200 kpc from the host (with random 
angular coordinates) and 0 and 200 kms“^, respectively. 
The model is randomly rotated and observed along the z 
direction — i.e. the z direction is taken to be the line-of- 
sight and the xy plane is taken to be the plane of the sky. 

(ii) Disc model — In this model, 2 — 5 satellites are placed 
randomly between 0 and 200 kpc from the origin on the xy 
plane (prior to any rotation taking place) and then randomly 
given a z coordinate between —10 and 10 kpc. All satellites 
are assigned a 3D velocity of 100 kms“^, such that each 
satellite is in circular motion about the host, initially rotat¬ 
ing in the xy plane. The model is then randomly rotated and 
viewed along the z axis. Finally, to mimic observational error 
in the line-of-sight velocities, we add to the z component of 
each satellite’s velocity a random offset drawn from a normal 
distribution with a standard deviation of av = 20 kms“^. 


^ In particular, the toy models do not account for potential ve¬ 
locity correlations due to group infall, which [Cautun et al.| ( |201^ 
argues may be important, nor do they capture the complexities 
of observing against a background of interloper galaxies with po¬ 
tentially corellated velocities. 

^ In each case, where the model is permitted to have more than 
two satellites, we set the probability of a host having n satellites 
to be four times greater than the probability of having n + 1 
satellites. 
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Our qualitative and quantitative results are not strongly de¬ 
pendent upon the assumed velocity structure or thickness of 
the model satellite discs; since we are only concerned with 
the sign (+ or -) of the z component (post-rotation) of the 
satellite’s velocity vector, the magnitude of that vector - and 
any radial dependencies it might have - is largely unimpor¬ 
tant. 


(iii) M31 model — This model is based on the position and 
velocities of the 13 satellites belonging to the co-rot at ing 
plane identified around M31 by Ibata et ah (2013). The 
three-dimensional positions of the satellites are taken from 


McConnachie (2012) and the line-of-sight velocities are com¬ 


piled from McConnachie (2012) and Collins et al. (2013). 
Note that we only consider the 13 satellites exhibiting co¬ 
herent rotation; the two satellites aligned with the planar 
structure, but with counter-aligned line-of-sight velocities, 
are excluded. We assign each mock satellite a velocity, such 
that the radial component of its velocity is consistent with 
the observed line-of-sight velocities of the true M31 satel¬ 
lites, such that each satellite’s total velocity puts it in circu¬ 
lar motion around the host We then randomly select 2 — 5 
satellites to mock observe (independent of the luminosity 
of the true M31 satellites), and the system is randomly ro¬ 
tated and viewed along the z axis. We again add to the ^ 
component of the velocity a random offset drawn from a 
normal distribution with ay = 30 kms“^, representative of 
measurement error in the line-of-sight velocity. While con¬ 
structed using the positions and velocities of the co-rotating 
satellites in the M31 plane, it is useful to note that our model 
is not a true analog of the observed system as every member 
of the M31 planar structure is fainter (by ~ 1 —2 magnitudes 
or more) than our satellite luminosity limit (Mr < —16). 


(iv) Dumbbell model — In this model, each host is restricted 
to exactly two satellites. Once the first satellite is randomly 
placed on the xy plane (again, prior to any rotation), the 
placement of the second satellite is restricted, such that the 
opening angle between the two satellites is less than 10° 
when the system is viewed along the z axis. From there, each 
satellite is assigned a z coordinate between —10 and 10 kpc 
and the system is subject to random rotation, assigned of 
line-of-sight velocity errors, and finally viewed along the z 
axis. In essence, this model requires two satellites to be on 
opposite sides of their hosts and orbiting in rigid-body rota¬ 
tion. As was the case for the disc model, our results do not 
strongly depend on the adopted thickness of the dumbbell. 


For each realization of a model, we rotate the system 
to a random orientation before observing along the z axis. 
In our analysis, we simulate a variety of statistical sam¬ 
ples, each consisting of iV = 10® model realizations, where 
most samples include a mix of realizations drawn from the 
isotropic model along with one of the other three models. 
For example, the 50% disc + 50% isotropic model consists 
of 5 X 10® realizations of the disc model and 5 x 10® realiza¬ 
tions of the isotropic model. Note that such a sample does 
not consist of 10® hosts whose satellites have a 50% prob¬ 
ability to be placed in a disc and a 50% probability to be 


The origin prior to rotation is taken to be the position of a 
Milky Way observer. 


Model 

Psat 


X^ 

P 

100% Disc 

1.0 

460.07 

28.76 

< 0.001 

50% Disc -h 50% Isotropic 

0.71 

179.77 

11.39 

< 0.001 

25% Disc + 75% Isotropic 

0.50 

71.20 

4.45 

< 0.001 

10% Disc -h 90% Isotropic 

0.32 

22.44 

1.40 

0.13 

50% M31 + 50% Isotropic 

0.71 

221.79 

13.05 

< 0.001 

50% Dumbbell -h 50% Isotropic 

0.71 

8.70 

0.58 

0.89 

10% Dumbbell + 90% Isotropic 

0.32 

21.62 

1.44 

0.12 

100% Isotropic 

0 

11.40 

0.67 

0.83 


Table 1. reduced y^, and p values based on a comparison 
of the observed fraction of co-rotating satellite pairs in the SDSS 
versus that for various statistical samples of model satellite dis¬ 
tributions as described in ^ Models in which a large fraction of 
hosts harbor discs of satellites (including the model based on the 
M31 plane) are disfavored relative to our dumbbell model or the 
simple isotropic case. Calculations are made taking satellite pairs 
over the full range of a (i.e. 0° < a < 180°), but the p-values 
are largely unchanged when we restrict the comparison to a nar¬ 
rower range of opening angles (e.g. a < 35°). Also shown is the 
probability that an individual satellite is found in a disc traced 
by bright satellites, Psat- 


placed randomly. While we do not explicitly explore cases of 
this type, they are equivalent to the cases we explore that 
have a percentage disc composition of approximately Ps^t? 
where psat is the probability of an individual satellite being 
placed in a disc. Note that this is only an approximation, 
as some hosts in our models have three or more satellites. 
For example, a case where satellites of each host indepen¬ 
dently have a 50% chance of being placed in a disc, would 
be equivalent to our 25% disc + 75% isotropic model. 

In Table 1, we present a summary of our analysis of the 
y^ values describing the goodness of fit for several statistical 
samples of modeled systems in comparison to the observed 
fraction of co-rotating satellite pairs as a function of opening 
angle in the SDSS (see Fig. If)- The number of degrees of 
freedom that enter into each p-value calculation is based 
on how restrictive the model under consideration is: disc 
models are taken to have one fewer degree of freedom than 
the purely isotropic model or the model that is based on 
the observed positions of M31 satellites, since satellites are 
restricted to being placed in the disc. The dumbbell model 
essentially introduces an additional constraint on the disc 
model, so dumbbell models have yet one fewer degree of 
freedom. 


3.1 Disc Model and M31 model 

In Figure]^ we show the observed fraction of SDSS satellite 
pairs that are co-rotating as a function of the opening angle 
a in comparison to mock observations of various disc models. 
The thick blue line corresponds to a statistical sample com¬ 
prised purely of satellite discs, while the green, orange, and 
magenta lines correspond to samples composed of 50%, 25%, 
and 10% satellite discs, respectively, with the remainder of 
each sample consisting of isotropic satellites. In addition, 
the dashed green line corresponds to a sample with satellites 
for 50% of the simulated hosts following a M31 model and 
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the other 50% distributed according to an isotropic satellite 
population. For comparison, the solid black line shows the 
co-rotating fraction for a purely isotropic sample. 

We find strong disagreement between the fraction of 
co-rotating satellite pairs as a function of opening angle for 
models in which 25% or more of hosts harbor planes of satel¬ 
lites in comparison to that for the observed SDSS sample. 
In particular, the presence of inclined planes (relative to the 
line-of-sight) in the toy models results in a significant over¬ 
abundance of co-rotating pairs at 20° ^ < 60° in com¬ 

parison to the SDSS observations. Overall, the SDSS data 
agree reasonably well with models where at most ~ 10% of 
hosts have satellites residing in planes (or ~ 90 — 100% of 
the hosts have satellites distributed isotropically in phase 
space); however, the 100% isotropic model does fail to re¬ 
produce the overabundance of co-rotating pairs at very small 
opening angles {a < 10°), as measured in the SDSS sample. 
As shown in Fig. the M31 sample follows the 50% disc 
+ 50% isotropic sample very closely, as they fundamentally 
represent the same satellite arrangement (with the caveat 
that the positions and velocities of the satellites in the M31 
model are tailored to match the observed positions of the 
true M31 satellites). While our results suggest that coher¬ 
ently rotating discs of bright satellites are not common, the 
objection could be raised that the velocity selection criteria 
used to select the SDSS systems systematically removes in¬ 
clined satellite planes (i.e. those systems with satellite pairs 
at opening angles of 10° < a < 170°); we address this pos¬ 
sible selection effect in §3.2| 

3.2 Velocity Modeling and Cuts 

As highlighted in §3.11 those toy models, in which a high 
fraction of host galaxies harbor satellite planes, are disfa¬ 
vored in part due to the lack of excess co-rotating satellite 
pairs at intermediate opening angles (i.e. 20° < o < 60°), 
corresponding to satellites in discs at non-zero inclination 
angles. If our sample selection criteria, in particular our ve¬ 
locity cut, are biasing us strongly against such systems, we 
could perhaps reconcile the apparent discrepancies between 
the disc models and the SDSS data. In a simplified test case, 
where each satellite orbits their host in a disc at a velocity 
of Vb, imposing a velocity cut of exactly Vo on the host- 
satellite velocity offset would retain only perfectly edge-on 
discs, leading to a signal much like the one observed at small 
opening angles. Relaxing this velocity cut would permit pro¬ 
gressively more face-on discs, and imposing no velocity cut 
would in principle permit any disc inclination angle. In con¬ 
structing our model, we assigned a characteristic velocity of 
100 kms“^ to the satellites and only selected those satellites 
that have a ID velocity offset (relative to their host) greater 
than X 25 km s“^. Since the toy model is essentially scale- 
free, this is equivalent to removing satellites that have a ID 
velocity offset less than \/2 x 25% of the characteristic 3D 
velocity for satellites of the host — i.e. a velocity threshold 
of 0.35 Vb is applied, where Vb is the characteristic 3D veloc¬ 
ity of the satellites. The higher this velocity cut, the more 
we would expect a contribution to the co-rotating fraction 
at intermediate a from inclined discs to be suppressed, as 
such systems would have a significant component of their 
velocity moving perpendicular to the line-of-sight. 

In Figure we demonstrate the impact of altering our 



Figure 5. The fraction of co-rotating satellite pairs in the SDSS 
satellite sample (cyan line) and in the 25% disc + 75% isotropic 
model (red lines), varying the value of the velocity selection limit 
applied to the models. The dotted, dashed, and solid red lines 
correspond to models in which only systems with AViog > 0-3, 
0.6, 0.9 Vo are included, respectively, and where Vq denotes the 3D 
characteristic velocity of satellites in the plane/disc. The dramatic 
overabundance of co-rotating pairs at small a in the SDSS data is 
only captured by models that contain nearly exclusively edge-on 
discs (i.e. applying the AVios > 0.9 Vq cut). 


velocity selection criterion on the measured co-rotation frac¬ 
tion versus opening angle for the 10% disc + 90% isotropic 
model. The various red lines range from selection limits of 
30% of the 3D velocity to 90% of the 3D velocity, with the 
SDSS data included for comparison (cyan line). Reproduc¬ 
ing the lack of excess co-rotation at 20° < < 60° — or 

similarly the sharpness of the increase in co-rotating satel¬ 
lite pairs at very small a — requires a very high velocity 
cut (i.e. ~ 0.9 Vo). Given that our adopted velocity limit is 
Aldos > \/2 X 25 kms“^, such a strong selection bias would 
require that the characteristic velocity (Vb) of satellites in 
planes would need to be ~ 40 kms“^. If the planar satellites 
have such low velocities, a velocity cut of \/2 x 25 kms“^ 
would indeed correspond to 0.9 times the characteristic ve¬ 
locity of the plane members, and we could confidently state 
that we had removed all but edge-on planes. 

The “toy model” adopted by 114 as a comparison 
to their measurements of the co-rotating satellite fraction 
in the SDSS utilizes exactly this characteristic velocity 
(40 kms“^); as a result, their toy model only includes edge- 
on (or nearly edge-one) discs, thereby reproducing the over¬ 
abundance of co-rotating satellite pairs at small opening an¬ 
gles (see their Fig. lb). Assuming a 3D characteristic veloc¬ 
ity of 40 kms“^ for satellites in planes, however, is inconsis¬ 
tent with the broad expectations from subhalo kinematics 
in iV-body simulations as well as the observed line-of-sight 
velocities of satellites in the M31 plane, which have a me¬ 
dian value of \AVios\ ~ 92 kms“^. Moreover, such a low 
characteristic velocity directly conflicts with the measured 
line-of-sight velocity offsets of the co-rotating satellite pairs 
identified by 114 (see their Table 1) as well as those in our 
sample. 

Figure shows the host-satellite (line-of-sight) veloc- 
© 2015 RAS, MNRAS QQQ,[HfTol 




Are planes of satellites ubiquitous? 7 



Figure 4. The fraction of co-rotating satellite pairs as a function of opening angle (a) for various versions of the disc model, including 
the model based on the M31 plane, in comparison to the observed co-rotating fraction measured in the SDSS (cyan line). The observed 
kinematics of bright satellites in the SDSS are consistent with at most 10% of hosts having disc-like or planar satellite distributions, with 
the remainder being distributed isotropically. 
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Figure 6. The line-of-sight velocity offset for each satellite in the 
SDSS sample relative to its host as a function of its projected dis¬ 
tance from the host (red circles). Those satellites corresponding 
to co-rotating pairs with opening angles of a < 10° are high¬ 
lighted as cyan squares. Given that most satellites in the sample 
(especially those that are co-rotating and at small opening an¬ 
gles) have velocity offsets greater than the lower selection limit, 
we conclude that our observational sample is not strongly biased 
against inclined planes (or discs). 


ity offset plotted against host-satellite projected distance 
for each satellite in our SDSS sample, highlighting those 
systems belonging to pairs with a < 10°. The mean ab¬ 
solute value of the ID (line-of-sight) velocities for the 44 
satellites in co-rotating pairs with a < 10° is 109.4 kms“^, 
consistent with that for the overall sample (111.4 kms“^). 


This indicates that a velocity cut of y/2 x 25 kms“^ would 
be insufficient to select only edge-on planes, such that we 
should also detect co-rotation from planes at moderate in¬ 
clination angles (i.e. yielding an elevated co-rotating fraction 
at 10° < a < 60°). As a result, we argue that the observa¬ 
tional data are inconsistent with co-rotating planes being 
ubiquitous in the local Universe — at least with respect to 
satellites brighter than Mr = —16. 


3.3 Dumbbell model 

Recognizing that the observed variation in co-rotating satel¬ 
lite fraction with opening angle is strongly inconsistent with 
discs or planes of satellites being prevalent, we now discuss a 
potential alternative scenario: the dumbbell model (see it- 
Figurej^shows the co-rotating satellite fraction as a function 
of opening angle for the SDSS satellite sample in comparison 
to two formulations of the dumbbell model, one with 50% 
dumbbells (and 50% isotropic satellites) and one with 10% 
dumbbells (and 90% isotropic satellites). Both models are 
in relatively good agreement with the observational data. In 
particular, the dumbbell model is able to reproduce the ob¬ 
served overabundance of co-rotating pairs at small opening 
angles and the corresponding sharp decrease at a ~ 10°. 
In Q we address the physical motivation (or lack thereof) 
for dumbbell-like satellite conhgurations. For now, we note 
that the dumbbell model yields a better ht to the observed 
kinematic data than the disc model. The p-value for the 
toy model with 10% contribution from dumbbells is 0.948, 
meaning we fail to reject it at any conhdence level. We also 
fail, at the 90% conhdence level, to reject the 50% dumbbell 
model {p = 0.152). While we similarly fail to reject the 10% 
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Figure 7. The fraction co-rotating satellite pairs as a function 
of opening angle a for the observed SDSS sample (cyan line) 
and various simulated dumbbell samples (green, magenta, and 
black lines). While unlikely to be physically accurate, a model in 
which a small fraction 10%) of massive galaxies host satellites 
in dumbbell-like distributions is able to reproduce the overabun¬ 
dance of co-rotating satellites pairs at small opening angles. 


disc model at 90% confidence {p = 0.223), we reject all other 
disc models at > 99% confidence (see Table |^. 


4 DISCUSSION 

Our analysis is motivated by the work of 114, which found 
an increase in the incidence of co-rotation in satellite pairs 
very near diametric opposition with respect to their host. If 
this excess co-rotation is the result of physical processes, it 
would tell us something significant about the co-evolution 
of satellite systems and the behavior of dark matter on 
small scales - perhaps indicating a serious problem with 
the ACDM model. On the other hand, we must consider the 
possibility that the co-rotation signal of 114 is not robust 
and is the product of random chance. In this section, we 
weight these competing possibilities. 

In selecting our sample, we adhere closely to the selec¬ 
tion criteria of 114; however, we did not reproduce the 114 
satellite sample exactly. Of the host systems listed in Table 1 
of 114, two galaxies fail our selection criteria. In one of these 
cases, the host magnitude fell just outside of our selection 
window, while the other host had a neighbor of nearly equal, 
but slightly brighter magnitude {AMr = 0.01). In addition, 
our selection includes several systems not identified by 114, 
though we do not expect these differences between our sam¬ 
ples to bias our results in any meaningful way. Of the 22 
satellite pairs with a < 8° in 114, 20 (91%) appear in our 
sample. Moreover, we reproduce the principal result from 
114, finding an excess of co-rotating satellite pairs located 
on diametrically-opposed sides of their host galaxy (i.e. at 
small opening angles). At a < 8°, we find 21 out of 27 satel¬ 
lite pairs to be co-rotating, corresponding to a co-rotating 
fraction of 0.78 zb 0.08%. 

However, when examining the kinematics of the entire 
population of satellite pairs (i.e. at all opening angles), the 


dependence of co-rotating fraction on a is inconsistent with 
the expectations from simple models of rotating discs (or 
planes) of satellites. For example, at large a, we would 
naively expect an overabundance of co-rotating satellites 
mirroring that detected at small opening angles. |Ibata et ah] 
(2014b) argue that caution must be taken when considering 
satellites on the same side of the host, as relative motion of 
the satellites with respect to each other (e.g. in orbital binary 
pairs) could dominate the motion of a satellite around the 
host, particularly when viewed in projection. As a result, 
when examining satellite pairs with large opening angles. 


Ibata et al. (2014b) require that the brighter of the satel¬ 


lite pair be at least two magnitudes fainter than the host 
to minimize the impact of infalling sub-groups; to mitigate 
deblending problems, identified pairs closer than 25 kpc in 
separation are also excluded. With these cuts in place, the 
114 sample contains 15 pairs of satellites with a > 172°, 
10 of which are co-rotating (Ibata et al. 2014b). Applying 
identical cuts to our sample, we find 13 out of 21 satellites 
co-rotating over the same range in a, while removing these 
additional selection criteria yields 44 pairs of satellites, 23 
of which are co-rotating. In each case, the satellite sample 
at large a is divided nearly evenly between co-rotating and 
counter-rotating pairs, and the co-rotating excess at large a 
is significantly less substantive than that at small a. This 
agrees well with the work of Cautun et al. (2015), which 
also finds no evidence of excess co-rotation in satellite pairs 
located on the same side of their host. 

The differences between the expected kinematics of 
satellites belonging to discs (or planes) and that observed 
in the SDSS satellite sample extend far beyond the large a 
regime. If the observed overabundance of co-rotating satel¬ 
lite pairs at small a is associated with a large population 
of satellite planes, then we would expect that planes at 
non-zero inclination angles to contribute co-rotating pairs 
at intermediate opening angles to the observed sample, re¬ 
sulting in a gradual decrease in the co-rotating fraction at 
0° < a < 90° (see Fig.[^. Instead, what is observed is a pre¬ 
cipitous drop in the co-rotating fraction, with 21 out of 27 
(78 zb 8.0%) satellite pairs found to be co-rotating at a < 8°, 
but only 27 out of 63 pairs (44 zb 6.2%) co-rotating over 
8° < a < 28°. Moreover, the dependence of co-rotating frac¬ 
tion on a is quite noisy (see Fig.|^, such that excesses in the 
co-rotating fraction - comparable to that found at a < 10° 
- are detected at other relatively narrow ranges of open¬ 
ing angle. For example, 20 out of the 31 satellite pairs with 
70° < a < 80° are co-rotating, which (if significant) is diffi¬ 
cult to reconcile with the predicted kinematics of satellites 
in planes. The discrepancies between the data and the disc 
models are borne out in the statistical analysis, with models 
involving more than 10% contribution from discs (or planes) 
strongly disfavored (see Table 1). This is in good agreement 
with the predictions of [Cautun et al. (2015), which argue 
that only ~ 15% of bright satellites should share a coher¬ 
ent sense of rotation to within 25° , as would be required for 
planar co-rotation. Our analysis strongly disfavors the claim 
that planar configurations of satellites are ubiquitous over 
the magnitude range {Mr < —16) considered. 

While the disc model does a relatively poor job of repli¬ 
cating the observational data, our toy model with satellites 
arranged in dumbbell configurations provides a better fit to 
the observed satellite kinematics. By construction, the satel- 
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lites in the dumbbell model are located at opposition, such 
that they able to replicate the sharpness of the increase in 
co-rotating fraction at small a while being consistent with 
no excess co-rotation at larger opening angles. The physi- 
cality of the dumbbell model is questionable, however, as 
it is difficult to imagine a scenario that involves satellites 
co-rotating in such a highly-constrained configuration. Ad¬ 
ditionally, the satellites belonging to M31’s plane are not 
in a dumbbell configuration. For these reason, we conclude 
that the data are not likely to be the result of dumbbell 
configurations either. 

In the absence of viable alternative models, we argue 
that the excess of co-rotating satellite pairs at small a is very 
likely the result of random noise. Among the disc models, 
only the model with a 10% contribution from satellite discs 
is a comparably good fit to the SDSS data as the purely 
isotropic model. Moreover, 25% of random realizations of an 
isotropic model comprised of 400 hosts (i.e. comparable in 
number to the observed SDSS sample) yield excess co- and 
counter-rotating pairs of satellites comparable in significance 
to that measured in the SDSS sample, although just ~ 1—2% 
show an overabundance of co-rotating pairs at small opening 
angles {a < 20°). 

These results should not be taken to say that there 
are no planes of satellites in the Universe. Statistically, the 
observational data are roughly consistent with < 10% of 
isolated, massive galaxies playing host to planes, or more 
specifically planes with multiple luminous satellites. Our re¬ 
sult does not exclude the possibility that planar structures 
preferentially populated by faint satellites could be ubiqui¬ 
tous; if true, this could provide powerful insight into the for¬ 
mation of planar satellite structures. With the caveat that 
all of the satellites under our consideration here are brighter 
than the satellites belonging to the M31 plane, we exclude 
the possibility of ubiquitous planes analogous to the M31 
plane at a very high confidence level. 


5 CONCLUSIONS 


In this work, we investigate the ubiquity of co-rotating 
planes of satellites similar to that observed around M31. 
Using data drawn from the SDSS, we study the orienta¬ 
tion and kinematics of bright (Mr < —16) satellite pairs 
around isolated galaxies, selected to be roughly analogous 
to the Milky Way and M31. By comparing the fraction of 
co-rotating satellites pairs as a function of opening angle to 
the predictions of simple models, we investigate the signa¬ 
tures of coherent rotation arising from satellites arranged in 
planar structures. Our findings are as follows: 


• We confirm an excess of co-rotating pairs of satellites at 
opening angles of a < 10° — i.e. satellite pairs located on 
diametrically-opposed sides of their host (see Fig. [^. This 
overabundance of co-rotating pairs at small opening angles, 
as first identified by 114, has been cited as evidence that 


50% of satellites belong to coherently-rotating planes (Ibata 
et al.||2014b| |. 


• We find that the excess of co-rotating pairs at small open¬ 
ing angles is unlikely to be due to ubiquitous co-rotating 
planes of satellites. While the behavior at small a is sug¬ 
gestive of planes aligned along the line-of-sight (i.e. with an 


orbital inclination of zero), the signal is strongly inconsis¬ 
tent with mock observations of satellite galaxy planes (or 
discs). In particular, we find no contribution to the co¬ 
rotating fraction from planes inclined relative to the line- 
of-sight (i.e. with satellites configured at opening angles of 
10° ^ < 60°). For our sample of isolated systems, we find 

that at most ~ 10% of hosts harbor satellites planes - as 
traced by the luminous (LMC-like) satellite population. 

• The excess of co-rotating pairs of satellites at small a is 
better fit by a “dumbbell” model, where satellites have co¬ 
rotating partners located opposite the host galaxy (but not 
a true plane). However, this model is very likely unphysical. 

• We do not rule out the possibility that ~ 10% of hosts 
having co-rotating planes of satellites, or correspondingly ~ 
30% of satellites residing in such planes. This case is similar 
enough to the isotropic case that the observed SDSS co¬ 
rotating fraction as a function of opening angle is roughly 
consistent with < 10% of massive host galaxies harboring 
planes of satellites. 
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